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I. INTRODUCTION 

It has long been observed that stationary flows over an obstacle which cross the speed of 
surface waves (Froude number F n = 1) are often associated with an undulation, i.e. a zero- 
frequency wave with a macroscopic amplitude (TJ [2] . In this paper we study a time dependent 
phenomenon which can be considered, in an appropriate zero-frequency limit, as a precursor 
of this phenomenon. We study the linear scattering of low frequency waves which propagate 
against a non turbulent and irrotational flow when the free surface does not already contain an 
undulation. When approaching the locus where F n = 1 is crossed, these waves are blocked while 
their wave vector increases. Using the analogy with gravity [3H7j, these flows correspond to 
white hole geometries, and the horizon is located where F n — 1. To understand the scattering 
of waves on such flows one must take into account the dispersion relation of gravity waves. The 
relation between their frequency Q (measured in the fluid frame) and their wave vector k is 

= c^Jk) = gkt a nh(h B k), (1) 

where fig the height of the background flow, and g the gravitational acceleration. In the 
usual case, when the conserved frequency u (measured in the lab frame) is high enough, an 
incoming wave packet engenders two reflected waves which follow different trajectories, and 
which possesses different relative amplitudes. In the limit u — >■ 0, two unusual properties 
combine and produce a long lasting undulation with a large amplitude. One first finds that the 
two scattered waves merge and form a standing wave with nodes at definite places. Second, one 
finds that the amplification of this wave grows as This second aspect is directly related 
to the famous prediction of Hawking [S] according to which an incipient black hole should emit 
a steady thermal spectrum of photons, see [9j [TO] for a pedagogical introduction and a review. 
In this sense black hole emission and the generation of undulations are based on a common 
amplification mechanism. 

To demonstrate that free surfaces (which contain no undulation) are unstable with respect 
to initial conditions, we consider a series of wave packets with decreasing value of u which 
all produce the same undulation. We show that the amplitude of the incoming wave linearly 
decreases in u. Hence there is an instability at very low frequency. So far this concerns the 
scattering of initial wave packets. However, given that in realistic conditions low frequency 
waves are inevitably excited in an uncontrollable manner, a statistical treatment of incoming 
wave amplitudes is required. When using a thermal ensemble to characterize their statistical 
properties, the above amplification is further enhanced since the noise in a thermal bath grows as 
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This extra power of 1 /lo leads to a time dependent growth of the undulation amplitude [TTI - 
IT3] . The growth will continue until some mechanism, which could be due to nonlinearties or 
dissipation, will saturate it. These late time processes are not discussed in the present work. 

It should be emphasized that similar phenomena occur in all dispersive media. In fact they 
are closely related to the appearance of the "layered structures" found in 4 He [H] , and in Bose 
gases [15], when the flow exceeds the Landau critical velocity 1 . They also occur in atomic Bose 
Einstein condensates where the dispersion relation is F^ EC = k 2 + £ 2 fc 4 /4 instead of Eq. Q. 
In this paper, even though we have in mind gravity waves, we formulate the problem in terms 
which apply to the general case. Namely, our treatment depends on the dispersion relation F 2 
which characterizes the medium, and quantities which characterize the background flow: the 
velocity of the flow, the speed of low wavenumber waves, and the fluid density. 



II. ACTION FORMALISM AND WAVE EQUATION 

We consider surface waves which propagate in a water tank of constant transverse dimension 
L±. We assume that the flow is incompressible, non turbulent, and irrotational. We also 
assume that both the bottom of the tank and the background free surface do not depend on 
the transverse coordinate, and become asymptotically flat in the upstream region. We call 
/i as and v as the asymptotic values of the water depth and the background flow velocity in this 
region. For simplicity, we only study waves with no dependence in the transverse coordinate 2 , 
and we neglect the effects of capillarity. To incorporate them, one should consider the dispersion 
relation which generalizes Eq. Q by including the capillary length [I]. 

When the background flow is stationary and non-uniform, the linear equation which governs 
the propagation of these waves is rather complicated (U [7]. Yet, this equation can be derived 
from an action for the perturbations <fi of the velocity potential. Moreover this action possesses a 
rather simple structure which is very similar to that of sound waves in an irrotational fluid [16] . 
As we shall see, the properties of the scattering of waves are efficiently described using the 
action formalism. The action has the following structure 

S=\f p{x) {^yP* + V M ^)0] 2 + 4>F 2 (d A (x)) J dxdt. (2) 

The function p(x) is an effective 1-dimensional fluid density, v x (x) is the background flow 
velocity, c(x) fixes the low frequency group velocity, and F 2 is a differential operator which 



1 We are grateful to T. Jacobson to have drawn our attention to these works. 

2 Modes with a non zero transverse momentum p± are studied in [13] . 
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governs the dispersion relation. The combination pF 2 forms a self-adjoint operator, and d\(x) 
is the local dispersive wavelength. (In an atomic Bose condensate, the latter is known as the 
healing length [T71 [18].) For each fluid, these functions are related in a specific manner to the 
properties of the background flow. 

For gravity waves, assuming an incompressible fluid, i.e., a constant 3-dimensional density 
Pq D , there are several (physically equivalent) ways to identify these functions. Using [7] 3 a 
convenient choice is 

F 2 (d A (x)) = 1 id x tanh(d A (x)id x ), (3a) 
u A {x) 

d A (x) = h as V ^^\ (3b) 
v z {x) 

p(x) = pl D L ± d A {x), (3c) 

c 2 (x) = d A (x) [g + d y (v 2 /2) y = hB{x) ] = d A (x) g cS (x). (3d) 

In the above, v 2 = v 2 + v 2 , where v x and v y are the horizontal and vertical components 
of the background velocity, evaluated along the free surface y = Iib(x). The quantity g c s is 
the effective gravitational acceleration which takes into account the centrifugal acceleration. 



Asymptotically, all x dependence are negligible, and Eq. (3d) delivers the standard expression 
C L = h &s g. For long wavelengths, i.e. low gradients kd A <C 1, the dispersive length d A drops 



out from Eq. (3a) and one gets the dispersionless expression F 2 — > k 2 , since k = —id x . For 
smaller wavelengths, combining F 2 and c 2 , one finds a generalized version of Eq. (jl]) where d A 
acts as a dressed value of h as . In fact, to lowest order in d y v x , d A {x) reduces to the water depth 
at x. 

As far as the scattering of waves is concerned, all is needed is the knowledge of the differential 
operator F 2 and the functions p, v, c, d A entering in Eq. (|2]). In other words, the intricate aspects 
of the above equations will play no role in the sequel. Yet, to make relevant physical predictions, 
one extra information is needed, namely, the relation between the velocity potential and some 
observable quantity. For gravity waves, the relevant observable is the vertical fluctuation of the 
surface 5h(t,x) = h(t,x) — h B (x) with respect to the background free surface y = h B (x). The 



link between 6 and 5h is 



5h(t,x) = l —-(d t + v x (x)d x )(j)(t,x). (4) 

g c ${x) 



3 In this work, Unruh derives the wave equation on stationary background flows using curvilinear coordinates. 
Here we re-express his results using horizontal and vertical cartesian coordinates x, y. The simplest way to 
obtain Eqs. p]), is to start from the wave equation in Eq. (113) and to use Eqs. (28), (29) to express in 
terms of d x . 
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It is interesting to notice that 5h is directly related to 7T, the momentum conjugated to <p, given 
the action of Eq. Indeed one has 

p(x) 



n(t,x) = ——(d t + v x (x)d x )(j)(t,x) 

C Z {X) 

= -(pl D L ± )6h(t jX ). 



(5a) 
(5b) 



It is also interesting to notice also that for sound waves, e.g. in an atomic Bose condensate, 
the density fluctuation 5p is related to the momentum 7T, and the potential <f), by very similar 
equations, see Eq. B.ll in [18J. Therefore the forthcoming analysis also applies to these waves, 
when using the appropriate F 2 operator governing dispersion. 

The conserved Hamiltonian associated with Eq. ^ will play an important role in the sequel. 
It is given by 

+ v 2 x (d x (f)) 2 - c 2 0F 2 (c/ A )0 } dx. (6) 



H 













v x d x <p 










P 



The wave equation can then be obtained from Hamilton equations. The first equation dt 



{0, H}, where { , } is the Poisson bracket, gives Eq. (5a). The second equation, d t 7i = {n, H}, 
gives 

{d t + d x v x (x))ir(t, x) = p(x)F 2 {d A {x)) </>(t, x), (7) 
which corresponds to Eq. (86) in [7j. Taken together, these equations give 



{dt + d x v x {x))4r\{dt + v x (x)d x ) - p(x)F 2 (d A (x)) 

C Z {X) 



(p(t,x) = 0. 



(8) 



Using the fact that the co-moving frequency Q is given by Q = u — vk, one verifies that the 
Hamilton- Jacobi equation associated with Eq. Q (obtained by writing = e - lU)t + l S k ^ dx ' an d 
by working to lowest order in the spatial gradients) gives 



n 2 = (oj - v x (x) k^x)) 2 = c 2 (x)F 2 (d A (x)). 
which generalizes Eq. ([T]) to non-uniform stationary flows. 



(9) 



A. Complete orthonormal mode basis 



When considering stationary background flows, as in the present case, it is extremely useful 
to decompose the surface waves in terms stationary complex modes of the form e~ tu)t ip^x) . The 
spatial part ^(x) then obeys 



(w + id x v x (x))^\(uj + iv x (x)d x ) + p(x)F(d A (x)) 
c z [x) 



(x) = 0. 



(10) 
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These modes are useful because they form a complete basis for the real solutions of Eq. ([8]), and 
because they efficiently describe their scattering. For more explanations, we refer to [18, 19J. 
The notions of completeness and orthogonality are based on the conserved product associated 
with Eq. Q. For any pair 0i, 2 of complex solutions of Eq. (J8|, the integral 



'* 1 7r 2 -ir* 1 (j) 2 )dx, (11) 



is constant, in virtue of Hamilton's equations. (This is true even when p, v , and c depend on 
time. When they are independent of time, the hamiltonian of Eq. ^ is another conserved 



quantity.) Several important properties concerning Eq. (11) should be mentioned. 

First, when considering two stationary modes, their overlap vanishes when uj\ ^ uj 2 . This 
guarantees that the orthogonal w-sectors do not mix with each other, and can thus be studied 
separately. 

Second, the sign of the above product is not definite. In fact the norm of the complex 
conjugated (e~^*0 w )* is the opposite of that of e~ ?w *0u>- To form a complete basis it is therefore 
sufficient to consider the set of positive norm modes together with their complex conjugated. 
Positive norm modes will be noted 0", and negative norm ones (<^)*, where a is a discrete 
index which accounts for the dimensionality (discussed below) of the set of modes at fixed u. 
As usual, it is appropriate to work with an orthonormal mode basis where all positive norm 
elements obey 



,) = N5 ab 8(u-u'), (12) 



where 5 ab is the Kronecker symbol, where 5 (u — u') is the Dirac distribution (because the 
domain of x is the entire real axis), and where N is an arbitrary (real and positive) constant 
which has the dimension of an action. A possible choice for gravity waves which depends on 
the flow properties is 

„3D j 

N gw = pl D L ± xhl s c as = P -^. (13) 

■ 1 i 1 — 3/2 

When using the scalar product (11) and the relation (5a), one verifies that iV = /i as Cas is the 



'net' amplitude of the modes (see Eq. (17) below), which guarantees that they are normalized 



as in (12). 



Third, the set of modes at fixed \u\, i.e. the set {</>", 6 _ tJ } and their complex conjugated, can 
be alternatively organized according to the sign of the frequency idt. This re-arrangement is 
justified because modes with opposite values of id t do not mix when propagating in a stationary 
flow. This means that e~ lult <t> a , will mix with e~ tult (<f) b _, ,)* but not with e tuJt d> b _, ,. Therefore, in the 
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sequel, we shall work with the set (0?L W )*} which contains modes with norm of both signs. 
In the presence of dispersion, and/or in inhomogeneous flows, there is no relation between the 
sign of u and that of the norm. Hence a case by case analysis is required to identify the various 
modes and their norm. 

Fourth, when considering the real solutions of Eq. ^ which vanish asymptotically (for 
x — > ±00), a complete basis is provided by the stationary modes W which are asymptotically 
bounded, i.e. stay finite as x — > ±00. Hereafter we call them ABM, for asymptotically bounded 
mode. Importantly, the dimensionality of this set depends on the value of u, and the properties 
of the flow. In what follows we consider flows which cross the speed of low frequency waves 
once, at a position conventionally put at x = 0, as described in Fig. [Tj In this case, there are 
three ABM for uj < w max , and only two for u > w max 4 The value of o; max is fixed by the condition 
that the line u — v^k be tangent to the dispersion relation evaluated on the asymptotic right 
side, see Fig. [2j Above a; max , we call the two orthogonal ABM by and (<jF_J)*, where the 
superscript U (V) indicates the sign of the group velocity with respect to the fluid: U (V) modes 
propagate against (with) the flow. Below w max , the three ABM are denoted <f>%, (<^ w )*, and 
4%, where the third one has a positive norm. Because we are interested in the low frequency 
limit, only the second set is relevant. 

To simplify the presentation, in the sequel, we shall neglect the l^-mode, and only consider 
the [/-modes 0^, (0_ w )* which are responsible for the undulation, and for the Hawking effect. 
However, as shown in [181 EQ], the coupling between 0^ and the [/-modes is non negligible. 
Yet, it does not affect the low frequency behavior of the mode mixing between 0^ and (0^ w )*: 
only some overall constants are affected. Since our purpose it to study the behavior of low 
frequency waves, we can, and shall, ignore the ^-modes. Under this approximation, dropping 
the superfluous index U, any low frequency upstream moving wave can be decomposed as 



■'max 



0(t, x) = 2Re I ^ e-»* [a w <j> w {x) + b u (<p- u (x))*] dwj , (14) 

where a u = (e~ lult 4>u\ <t>)/N and b w = — (e _:(W *(0_ w )*| <p)/N are the constant components of in 
the mode basis one has adopted, and where the minus sign in the last equation arises from the 
fact that (0-oj)* is a negative norm mode. 



Since the norm of Eq. (11) is identically zero for all real waves, it seems to be of small 



When the flow is everywhere subcritical F n < 1, for the dispersion relation of Eq. ([T]), at fixed u>, there exist 
four ABM, two with positive norm, and two with a negative one, as can be seen from Fig. [2j When the flow 
everywhere supercritical, F n > 1, there are two ABM with opposite norm. 



s 




Figure 1: Background profile function of x with v > 0. The chosen profile is v x — c = 

—cqD tanh (kx/cqD), where cq is the speed at x = 0. The frequency k is defined by d x (v x — c)\ x=0 = —k. 
In a gravitational context, this quantity corresponds to the surface gravity of the horizon [3]. The 
parameter D governs the extension of the near wave blocking region where c — v x ~ kx. On the 
asymptotic subcritical right side, v x (x) and c(x) reach the constant values v as and c as - 

relevance for surface waves 5 . Similarly the sign of the frequency seems not well defined since 
both signs of id t appear in real waves, both for the waves made with a w w or b u (<f>- u )*. To 
understand that there is a physically relevant, and important, information, let us consider the 



energy function. Using the mode normalization of Eq. (12), the energy transported by Eq. (14) 
is given by the value of Eq. ^ evaluated with = 0: 

1 fw miB 
E = H[4>]= -(4>\id t 4>) = N u [K| 2 - \bj 2 } doo. (15) 



We first verify that the choice of N used in Eq. (12) drops out from the products iV|a w | codoj 



and —N\b u \ 2 uj(kj, which have the dimension of an energy. We also see that these products 
give the energy contribution of real waves described by W , and that of those described by 



5 The norm of Eq. (Ill (with N = H) is used in quantum field theory to decompose the field operator into 
creation and destruction operators which provide a particle interpretation to the field excitations |21l 122) . In 
the case of phonons in a Bose condensate, this is well defined and also physically relevant [IT] . 
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respectively. (Hence working with the positive norm complex modes to describe 

the real waves, the sign of the frequency of the complex modes fixes the sign of the energy 
of the corresponding waves. Equivalently, when working with the positive frequency mode set 
the sign of the mode norm fixes that of the wave energy.) 



B. 5-matrix description of waves scattering 

The most efficient way to describe the scattering of low frequency waves in a stationary 
flow is to work at fixed u, and to study how the two complex modes W and (<j>- u )* mix. As 
in any scattering, it is useful to introduce the two incoming modes c/>™, (t/^)* which shall be 
scattered, and to relate them to the outgoing modes 0° ut , (0°^)* which have been scattered. 
Given that the flow is asymptotically uniform, the asymptotic properties of these four modes 
are easily found. Indeed, since the background functions v,c,p are all asymptotically constant, 



the solutions of Eq. (10) are superpositions of plane waves. Each plane wave is characterized by 



a constant wave number k w which is a real root of Eq. ^ evaluated asymptotically. The various 
roots are plotted on Fig|2j As explained earlier, we only consider those corresponding to waves 




Figure 2: Dispersion relation, f2 versus k, on the asymptotic right side and for the counter-propagating 
waves with respect to the fluid. The continuous (dotted) straight line is u — v as k of Eq. ^ evaluated 
asymptotically for to > (—to < 0) while the curve is c as F(k). Each intersection is thus a solution. 
The identification of three roots is described in the text. By comparing the straight and dotted lines, 
we see that k°^ > solves Eq. ^ for —to < 0, which means that — solves Eq. Q for to > 0. 



When considering the limit co — > 0, k° and merge and define kz of Eq. (16a). 
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moving upstream (with respect to the fluid, not the laboratory frame). On the asymptotic 
right side, where < w as < c as , they are three real roots. The low wave- vector one describes 
the incoming mode, whereas the two high wave-vector roots describe the outgoing modes. For 
uj — > 0, these roots reach equal and opposite values ±kz- Indeed to first order in uj one has 

C* ~ k z + u/v%, (16a) 
K ~ --4— = ^/<, (16b) 

where v g > is the group velocity v g = (d^k)' 1 evaluated asymptotically and for oj — > 0. In the 
second line one verifies that the in mode is propagating to the left since v g n = — c as + v as < 0. 
On the other side of the blocking region, for a; — > — oo where v > c, there is only one real 
root, which corresponds to the negative root — k^ u . To each root, there exists an asymptotic 



solution of the wave equation which can be normalized using Eq. (11). Of crucial importance, 
the norm sign of these modes is not the same. In fact, as we see from Eqs. (5a) and (11), 
this sign coincides with that of the co-moving frequency Q = u — vk. Hence, for u > 0, the 
two positive roots fc™ and /c° ut describe the in and out positive norm modes, whereas —k™ w 



and — k°^ describe the in and out negative norm modes. Explicitly, using iV of Eq. (13), the 
asymptotic positive and negative norm modes are 

. . / . Q. . . . 

fl = N 1 (17a) 

i^-S = N (17b) 
y/4ir\n(kZ)v%\ 

where the exponent a stands for in and out. (For (<p_ u )*, h as and c as should be evaluated on 
the left side when computing k 1 ^ and t>™, since this mode comes from that side.) 

We now define the globally defined in and out modes using the standard procedure [TU] . e.g., 



the positive ingoing mode (f)™ is the solution of Eq. ( 10 ) which asymptotes in the past (in the 



sense of a broad wave packet) to y?™. We use the symbol if to designate the asymptotic plane 



waves of Eq. (17), and <p the corresponding globally defined solutions. Since the two in modes, 



and the two out modes, form a basis of the ABM solutions of Eq. (10), they are related by a 
linear transformation which defines the S'-matrix 

\ / /A out \ 

(18) 

The fact that in each basis the two modes have opposite norm implies that the coefficients 
obey the anomalous scattering relation |aj 2 — 1/3^ | 2 = (ct^l 2 — 2 = 1 . Since \olJ\ > 1, the 
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scattering thus leads to an amplification of the waves [HI El 1231 [21] • In addition, in the limit 
uj — > 0, it can be shown that the amplification becomes very large, that is 

K| ~ > 1. (19) 

As we shall see, this is a crucial ingredient for the undulation. In [20] it was shown that, 
irrespectively of the smoothness of the profile near Froude number F n = 1, the coefficients 
diverge as 

~ \P„\ 2 ~ t£~- (20) 
When the profile is smooth, i.e. when nd^/co <C 1, one finds the Hawking result, i.e. R = k 
where k is the surface gravity defined in the caption of Fig. [TJ When the profile is so abrupt 
that the inequality is not satisfied, R is smaller than k and determined by the dispersive length 
g?a- To characterize the phase of the coefficients we write them as 

a.~K|e^ +9 ') and U ~ K|e^'-*\ (21) 

such that 9 is half the phase of and 9' a common phase. As we shall see, 9 participates 

to the determination of the location of the undulation nodes. 



III. SCATTERING OF WAVE PACKETS 
A. Incoming waves 

We consider a series of incoming wave packets of positive energy, sent from the right side 
against the flow. Using the in modes basis 0™, (0™ w )*, this means that we set = in 



Eq. (14). To get explicit expressions, we work with 



A ( (co-u)) 2 \ 

ex P - 0-2.-.2 > ( 22 ) 



raj V 2cT o w 

where A is a complex dimensionless amplitude A = \A\e l6 , and where the packets are normalized 
by 

la^du = \A\ 2 . (23) 



o 



We also assume that the waves are almost monochromatic. Irrespectively of the value of u, 
this is realized if 

Ob « 1. (24) 



In this regime, one finds n& = aoum 1 / 2 . The absence of any w-dependent phase in Eq. (22) 
guarantees that the packets are all centered around x = at t = 0. 
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We now analyze the early properties of this series of waves. Using the asymptotic behavior 



of 0™ given in Eq. (17) we get 



(x, t -oo) = 2Re { NA 1 - ^ 



e o" 



A/4vrwc as n y 2 



(25) 



In the broad wave regime of Eq. (24), we can accurately evaluate the integral with a saddle 



point approximation. Then, using Eq. Q asymptotically, where g c s = g, the corresponding 
incoming height variation is 



2-2 

5h in (x) ~ 6h% x sin (cot - k£x + 5) e -^ (t - x/v ^ 2 , 



where the amplitude is 



Sh 1 ; 



\A\N- 



2a 



(26) 



(27) 



As expected, the incoming wave (26) oscillates at a frequency u, has a wavenumber fc!?, and its 



envelope is propagating toward the wave blocking region at a speed v™ < 0. We see that the 
argument of A, the phase S, governs the precise initial positions of the nodes. We also see that 
at fixed A, its amplitude 5h m decreases linearly for with u. In this respect, it is also instructive 



to evaluate the conserved energy transported by the wave packet of Eq. (25). Using Eq. (15), 
one finds 

E = 2ooN\A\ 2 . (28) 
At fixed A, this energy is also arbitrary small in the limit uj — > 0. 



B. Outgoing waves 



At a time near t = 0, the waves of Eq. (25) reach the wave blocking region, around x = 0. 



There, they undergo a nontrivial scattering which is completely governed by Eq. ( 18 ) . Because 



of this mode mixing, two outgoing wave packets are generated and propagate to the right. 
To analyze these two waves, it is convenient to introduce the complex wave 0c, defined by 
6 = 2Re( 



in Eq. (14). At late time, using the out mode basis and Eq. (18), one finds 
4> c {x,t ->■ +oo) = 4>c(t,x) + <j>7,(t,x). 



(29) 



where the two waves are 
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where f2 out = Q(k° nt ) is the co-moving frequency of the outgoing modes, and t>° ut their group 
velocity. Evaluated through a saddle point method, one obtains 



2-2 

{x,t->+oo) ~ ^«^^° ut e-^e-^ ( *-^° ut ) 2 , (31a) 



0-( M _^ +oo ) „ ^^(^♦e-^e-VC*-^) 3 , (31b) 
where their common amplitude factor is given by 



A Q = AV2a u!7i 1 / 2 . (32) 
We now have all the ingredients to consider the limit u — > 0, keeping ctq constant. Using 



Eq. (21 ), the two outgoing waves merge to give a (real) zero- frequency wave of fixed profile, i.e. 



a^Te- lGlt + ^^%Te- iQt ~ 2\a Gl \e^ +e '^ u (x) } (33) 

where 

$ c/ (x) = Re{e 4 Vo Ut (^)}, (34) 

gives the profile of the undulation. Since asymptotically 0Q Ut (x) is given by Eq. (17), we get 



*u(x) = N C0 ^ X + d \ (35) 



We see that this profile is independent of t, u, and S. Using Eq. (33), the outgoing real wave 
neatly factorizes and takes the simple form 

2-2 

$(x,t +oo) ~ 4 |i4 fi >a ffi |cos((J + 9') x ^ u {x)e- Z ^-^ t - x ' v ^ 2 . (36) 

Since the slow varying envelope becomes a constant when u — » 0, the outgoing signal becomes 
a standing wave with nodes at fixed locations. This is our first result. 

To discuss the amplitude in physical terms, we compute the corresponding fluctuation of the 
free surface. Using Eq. (|4]), the undulation profile of Shu associated with Eq. (35) is 



Shu(x) = - V -^d x $u(x) ~ K S N, / V £ k ? sin (k z x + 9) . (37) 



Similarly, the surface wave associated with the packet of Eq. (36) is 

5h ont (x) ~ Sh ont sin (A; z x + 9) e - 5 % L (*- a, / , tf) 3 , ( 38 ) 

where the amplitude is 



5/i out = /i as iV J^r- \A Q a Q \ cos (5 + 9') . (39) 
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To get rid of the dependence on the amplitude \A\ and the normalization N, we compute the 



amplification factor. Using Eqs. (27) and (39), we find 



5h ont (c as -v as ) I v^kzK 

cos (5 + 9) . (40) 



5h in lo \/ 



Irrespectively of any choice, we see that the ratio diverges as 1/uj for u — > 0. This is our second 



result. We also see that the amplitude of the undulation critically depends on the cosine in (40) 
whose value is fixed by 5. In practice, this phase is hard to control with precision. We learn 
here that the resulting amplitude is drastically affected by a small change of initial conditions. 
Hence, unless one knows the initial conditions with an arbitrary precision, one is not able to 
predict what will be the amplitude of the undulation. The unpredictable character of the 
undulation amplitude (in the linear regime) is our third result. This property becomes even 
clearer when the state of incoming waves is a statistical ensemble. 

Before studying this, we make three extra comments which lead to specific predictions which 
could hopefully be tested in future experiments. 

First, we emphasize that the linear treatment we used predicts that both signs of the am- 



plitude are equally possible, since the sign is governed by the cos(5 + 9') of Eq. (40). 

Second, to be more specific, we suppose that t> as <C c as - This means that we work with deep 
water waves. It simplifies the expressions and is relevant for many experiments. In this case, 



using Eq. (16a), we have 

c 2 



kz = ^Y and v z = vj2. (41) 



Hence the net amplification factor of Eq. (40) becomes 



cos (5 + 9') . (42) 



We see that the amplification grows as for t> as — > 0. This growth will be regulated by 

capillary effects. 



Third, the phase 9 which governs the location of the asymptotic nodes in Eq. (37) is com- 
plicated to compute because it accounts for the propagation from the blocking region to the 
asymptotic region. On the contrary, the undulation profile has a rather universal behavior near 
the blocking region, where v x — c ~ — kx, that is, for \kx\ D (see Fig[T]). In addition we 
restrict our attention to the region where < l/d\ (i.e., \kx\ < 1), meaning that one can 
approximate F 2 {k) in Eq. ^ by k 2 — dj^k 4 /3. In this region, using the results of [12] and 
Eq. Q, up to an overall constant factor, we get 

5hu(x) oc Ai (-x/d hroad ) , (43) 
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4^=,^'] . (44) 



where Ai is the standard Airy function [25J, and where the effective length is given by 

6 k 

This length depends both the on dispersive length g?a(0) ~ /i.b(0), the height of the flow near 
the blocking point, and the gradient k = d x (v x — c) evaluated at that place 6 . In agreement with 
the results of [121 123 [28] , this broadening length governs the behavior of the undulation near 



the blocking point. We also note that the velocity potential $(/ of Eq. (34) is not described 
by an Airy function, but by its primitive integral. We finally note that these results also 
apply to undulations (of density fluctuation) found in Bose condensates where the dispersion 
is anomalous: -FJec = k 2 + £ 2 fc 4 /4. In that case, d\(0) is related to the healing length £ by 
<^a(0) = \/3£/2, and the undulation lives where the flow is supercritical. 

IV. STOCHASTIC ENSEMBLE 

In realistic conditions, low frequency modes are excited in a non controllable manner. To 
describe in simple terms this randomness, we assume that the incoming waves are character- 
ized by a thermal ensemble, i.e. with a noise growing like 1/u for oj — > 0. To describe this 



distribution, one must use the field decomposition of Eq. ( 14 ) in terms of in modes where the 
coefficients a w and b^ are random variables [29J. One sees again the usefulness of the mode 
basis. Such a Gaussian ensemble is defined by 

(a u ) = 0, (a*,a w ) = n"5(u - u/), (45a) 
(U = 0, (VM=nl6(u-«/), (45b) 

where the two variables are independent, i.e. {b^aj) =0. In a thermal state at temperature 
T, each monochromatic wave behaves as a (classical) harmonic oscillator in a thermal bath. 
However, since positive and negative energy modes respectively come from the right (R) and 
left (L) asymptotic regions, they will not share the same temperature. Moreover, their co- 
moving frequency Q = u — v k is Doppler shifted by the right or left asymptotic values of the 
flow velocity, see [IB] for a discussion about thermal states in fluid flows. To take both effects 

6 The evaluation of dbroad is most easily done by considering the dispersion relation ^ with F 2 ~ k 2 — d\k 4 /3 
and c One obtains nx ~ cok 2 d\/6, then using x = idk, the mode in fc-space is exp(id^ road fc 3 /3). 

Using a similar reasoning, one obtains the (w-dependent) broadening length of for deep water waves with 
non vanishing frequency u. 
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into account, we parameterize the low frequency powers as 

k B T R 



n, 



Nuj ' 
k B T L 



(46a) 
(46b) 



In this state, the average value of the surface perturbation vanishes, i.e. 



(5h(t,x)) = 0. 



(47) 



This is easily understandable using the wave packet analysis of the preceding section. Indeed, 



in a thermal state, initial phases (i.e. 5 in Eq. (26) for the incoming waves) are random. Hence 
when averaging over 5 in Eq. (38), the mean value vanishes. On the other hand, ((5h(t,x)) 2 ), 



the spread of Sh is non trivial. To compute it, we use the mode decomposition (14) together 



with the distribution (45 ) for the coefficients a u , b u . By a calculation similar to that of Eq. (33 ), 
and using the profile Shu of Eq. (37), one finds 7 



'nl + n b u )\aj\ 2 djj x (Sh v (x)f 



(48) 



This expression shows that the profile of the undulation is fixed and given in Eq. (37). In 



other words, the relative amplitude and the position of the nodes are not subjected to any 
randomness. On the other hand, its amplitude is a stochastic variable whose spread is fixed 
by the above expression. Therefore, the linearized treatment predicts that there is a high 
probability of observing a macroscopic undulation, with equal probability to find either sign. 



(Using Eq. (46), we verify that the choice of iV in Eq. (13) enters in Shu of Eq. (37) so that 



the spread of Eq. (48) is independent of N.) 



In addition, when integrating over low frequencies, the integral which fixes the spread of the 
amplitude diverges since its integrand behaves as 1/u 2 . To regulate it, we consider a flow which 
has been formed for a finite amount of time t. Since it is impossible to distinguish frequencies 



below ~ 1/t, we introduce a low frequency cut-off on the integral in Eq. (48) and consider 



k B {T R + T L ) 
2nN 



Kt. 



(49) 



We see that the diverging character for low frequencies engenders a linear growth in time. (This 
result has been confirmed by numerical simulations in atomic Bose condensates [U].) We also 
see that the low frequency waves coming from the right and the left both contribute, with 



7 For more details, we refer to the Chapter 4 of [3D], which presents results from [T21 [15] , 
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their respective powers given in Eq. (46). Of course, the amplitude of the undulation will grow 
linearly in time only as long as one can neglect nonlinearities, dissipation, or an infrared cut-off 
as in the case of transverse modes [13J. Ultimately, these will alter the picture and saturate the 
growth of the amplitude. 

V. CONCLUSIONS 

In this paper we studied the scattering of low frequency waves in supercritical flows when 
the stationary free surface is flat. In the zero-frequency limit, we showed that the outgoing 
reflected waves possess very specific properties. First, the two outgoing waves merge and form 



a single wave with a fixed spatial profile and nodes at specific places, see Eq. (33). Second, 



the amplification factor relating the amplitude of the incoming and outgoing waves diverges 



as the inverse of the conserved frequency, see Eq. (40). Third, this factor also depends on a 



combination of phases which renders the outcome effectively unpredictable, see the cos(5 + 9'), 



where 5 is given in Eq. (26) and 9' in Eq. (21). Fourth, we show that near the blocking point 



the profile of the undulation is given by an Airy function governed by a composite length scale 



formed with the dispersive length and the gradient of the flow k, see Eq. (44). 

These properties tell us that free surfaces (which contain no undulation) are unstable when 
sending low frequency incoming waves, and that the large and unpredictable amplitude of the 
undulation is an expression of this instability. To account for this unpredictable character, in 
the last part of this article, we computed the averaged properties of undulations when the low 
frequency incoming waves are described by a statistical ensemble. In agreement with [TTHT3"] . 



we found in Eq. (48) that the undulation spatial profile is not subject to any randomness while 
its amplitude is a random quantity. In addition, the spread of this amplitude diverges for very 
low frequencies. This divergence is regulated when considering that the stationary flow only 
existed for a finite lapse of time. The linear treatment predicts a growth which is linear in this 
lapse. It would be interesting to experimentally test this prediction, as that concerning the 



profile of the undulation given in Eq. (43). 

The low frequency divergence should also be regulated when including nonlinear effects. In 
this respect it would be particularly interesting to understand the transition from a growing 
undulation whose amplitude is a random variable to a saturated amplitude which is determin- 
istically fixed by nonlinear equations. This question is relevant for both classical fluids and 
quantum ones, such as dilute atomic Bose gazes. 
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